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Observed variations in the H ll region luminosity function (Hll LF) seen in spiral 

arm vs. interarm regions, and different galactic Hubble type, can be explained simply 

by evolutionary effects and maximum number of ionizing stars per cluster. We present 

Monte Carlo simulations of the Hll LF, drawing the number of ionizing stars A^^, from 

t^ ' a power-law distribution of constant slope, and the stellar masses from a Salpeter IMF 

•^ ■ with an upper-mass limit of 100 M©. We investigate the evolution of the Hll LF, as 

determined by stellar main-sequence lifetimes and ionizing luminosities, for a single 

burst case and continuous creation of the nebular population. Shallower H ll LF slopes 

measured for the arms of spiral galaxies can be explained as a composite slope, expected 

for a zero-age burst population, whereas the interarm regions tend to be dominated by 

evolved rich clusters described by a single, steeper slope. Steeper slopes in earlier-type 

galaxies can be explained simply by a lower maximum N^, cutoff found for the parent 
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;_i ■ OB associations. The form of the Hll LF can reveal features of the most recent {^ 10 

c/3 . Myr) star formation history in nearby galaxies. 
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1. Introduction 

H II regions have long served as a primary indicator 
of high-mass star formation, both in the Milky Way 
and in external galaxies. The luminosity function of 
H II regions (H II LF) is therefore a vital probe of the 
present global star formation properties in individual 
galaxies. In recent years, the number and accuracy of 
compiled HII LFs for nearby galaxies has increased 
substantially, extending the use of this tool to greater 
distances. 

The differential H II LF is usually parameterized as 
a power law: 



N{L) dL = A L-" dL 



(1) 



where N{L) dL is the number of nebulae with lumi- 
nosities in the range L to L + dL. Some interesting 
patterns have emerged regarding the form of the H II 
LF: 

1. A large number of galaxies show a break in 
slope of the HII LF with the fainter H II re- 
gions showing a shallower slope compared to the 
high-luminosity objects (e.g., Kennicutt, Edgar, 
& Hodge 1989, hereafter KEH; Rand 1992; 
Walterbos & Braun 1992; Rozas, Beckman, & 
Knapen 1996a;) 

2. The HII LF for arm and interarm regions of 
spiral galaxies sometimes show steeper slopes in 
the interarm regions (e.g., KEH; Banfi et al. 1993; 
Rand 1992). However, there are also many in- 
stances where no difference is found (e.g., Rozas 
etal. 1996a; Knapen etal. 1993; Knapen 1997) 

3. The slope of the H II LF is correlated with galac- 
tic Hubble type, such that early-type galaxies 
show steeper slopes than late-types. The power- 
law index a ^ 2.0 for Sb - Sc galaxies, compared 
to a ^ 1.7 for Sc - Im galaxies (KEH; Banfi 
etal. 1993). It is even steeper in Sa galaxies, 
with a - 2.6 (CaldweU etal. 1991). 

These differences in the HII LF have usually been 
interpreted as resulting from corresponding differ- 
ences in star formation properties. For example, the 
break in slope (Point 1, above) has been suggested 
to be caused by a physical transition between normal 
H II regions and the class of supergiant H II regions 
such as 30 Doradus in the Large Magellanic Cloud 
(KEH). Likewise, the difference in slope for arm vs. 



interarm regions, and among various Hubble types 
(Points 2 and 3), have been attributed to differences 
in gasdynamics and molecular cloud mass spectrum 
(e.g., Thronson, Rubin, and Ksir 1991; Rand 1992; 
KEH). However, we feel that it may be premature to 
draw a direct connection between the form of the H II 
LF and details of star forming environment such as 
these. There are other effects that can influence the 
form of the HII LF, that have not been adequately 
explored. 

One important effect, in particular, is the influence 
of simple evolution in the ionizing clusters, and hence, 
luminosity of the host H II regions. There are several 
empirical factors that suggest that nebular luminosity 
evolution is an important effect in the form of the H II 
LF. The steeper slopes seen in interarm regions of spi- 
ral galaxies (Point 2, above) have been identified by 
von Hippel & Bothun (1990) as evidence of an aging 
effect in those populations. Assuming star formation 
takes place primarily in the spiral density waves de- 
lineated by the arm regions, then the interarm regions 
should exhibit a typically older population of nebulae 
in the wake of this star formation activity. In addi- 
tion, von Hippel & Bothun cite the behavior of Ha 
equivalent widths as further compelling evidence of 
the importance of luminosity evolution. As a clus- 
ter ages, its associated Ha emission decreases, while 
that of the underlying red stellar continuum increases, 
thereby reducing the Ha equivalent width as a func- 
tion of time. They demonstrate that for the H II re- 
gions in NGC 628, there is indeed a clear decrease in 
Ha equivalent width with decreasing Ha luminosity, 
indicating higher numbers of evolved nebulae among 
the fainter populations. Finally, Knapen etal. (1993) 
and Rozas, Knapen, & Beckman (1996b) show that 
the slope of the nebular luminosity vs. volume rela- 
tion is slightly flatter in the interarm regions of five 
galaxies as compared to the arm regions. As empha- 
sized by the authors, this may be caused by obser- 
vational selection effects, but if real, we suggest that 
it is consistent with a higher tendency for shell-like 
geometries in the interarm nebulae, as would be ex- 
pected for an older population. 

All of these arguments suggest that evolved neb- 
ulae will be an important population in the H II LF. 
Von Hippel & Bothun (1990) carried out the first de- 
tailed investigation of the evolutionary effect on the 
HII LF. They constructed two models to reproduce 
the observed slope for NGC 628: one with zero evolu- 
tion, in which the power-law slope of the luminosity 



function is due entirely to the cluster mass spectrum; 
and one in which the distribution in nebular luminos- 
ity is due entirely to a different age and number dis- 
tribution of objects with the same initial cluster mass. 
The H II LF resulting from the evolutionary model is 
sensitive to the stellar initial mass function (IMF), 
since the relative numbers of stars at a given ion- 
izing luminosity become important when these stars 
become the dominant contributor to the nebular lu- 
minosity. Recognizing this, von Hippel & Bothun ori- 
ented their study toward exploring the use of the H II 
LF to place constraints on the IMF. As a result, they 
did not explicitly examine the consequences of ag- 
ing on the form of the HII LF itself. In addition, 
their assumption of uniform initial cluster masses in 
the evolutionary model is of only limited use, since 
actual initial cluster masses apparently have a power- 
law distribution (e.g., Oey & Clarke 1997; see below). 
Therefore, we will here investigate the effect of cluster 
aging on the form of the H II LF itself. 

2. Saturated Population 

We have carried out a preliminary investigation 
of this problem as part of a different study (Oey & 
Clarke 1997), and review the important relevant is- 
sues here. A first critical point is the fact that not all 
clusters have the same relative contribution to ioniz- 
ing luminosity from the various stellar masses. Such a 
situation is true for rich clusters that sample the IMF 
well, but breaks down for sparse clusters in which 
small-number statistics determine the relative ioniz- 
ing contribution of different stellar masses. We shall 
refer to the rich clusters with good stellar statistics as 
"saturated" with respect to the IMF, and those with 
poor statistics as "unsaturated." The importance of 
this fundamental difference in the dominant ioniz- 
ing contributors between the saturated and unsatu- 
rated clusters is demonstrated by a break in the pre- 
dicted slope of the unevolved HII LF. This change in 
slope was clearly demonstrated by McKee & Williams 
(1997; hereafter MW97), who performed Monte Carlo 
simulations of clusters drawn from a truncated power- 
law distribution in the numbers of member stars N^ , 
and a Scab (1986) IMF. Since L ex TV, for the satu- 
rated objects, the H II LF slope is the same as the par- 
ent slope in N^,] but for unsaturated populations, the 
HII LF slope is distinctly flatter. This flattening is 
caused by the increased scatter in the Us contributed 
by a bin of given A^* . The unsaturated condition oc- 
curs when this scatter is large compared to the mean 



L at the given A*, hence the slope of the HII LF 
flattens. MW97 observed that in their simulations, 
this change in slope occurs at Lsat = ^upj the lumi- 
nosity contributed by a single star at the upper-mass 



limit 



'■up I 



of the IMF. Stellar models currently sug- 



gest that ^up is in the range 38.0 < log^up < 38.5 
(Schaerer & de Koter 1997; Vacca, Garmany, & Shull 
1996; Panagia 1973). We caution however, that the 
slope break Lgat in the HII LF does not necessarily 
occur at ^up, as we show in §4.3. 

Investigation of evolutionary effects on the H II LF 
should therefore consider the differences in stellar ion- 
izing populations between the saturated and unsatu- 
rated regime. We examined the behavior of the sat- 
urated case analytically in Oey & Clarke (1997). We 
considered two extremes in the formation of ioniz- 
ing clusters: a single-burst scenario, in which all ob- 
jects are created at the same time; and a continuous- 
creation scenario. In what follows, we will also assume 
a constant IMF throughout. 

For the case of the single burst, it is apparent that, 
since all objects follow the same luminosity evolution, 
the initial slope of the HII LF is simply preserved. 

However, in the case of continuous creation, such a 
conclusion is not obvious. We described the luminos- 
ity evolution of the H II regions in terms of a fading 
function f{t) (Oey & Clarke 1997): 



L = Lofit) , 



(2) 



where Lq is the zero-age nebular luminosity. Stellar 
population synthesis work (e.g., Leitherer & Heckman 
1995; Beltrametti, Tenorio-Tagle, & Yorke 1982) sug- 
gests that the tail of the fading function f{t) is well- 
described by a power-law in time t, so that 
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t < t„ 



i/ims(mup) 



(3) 



t > ims(TOup) 



For saturated objects, the onset time, ims("T-up), for 
the fading is the same for all objects. Nebular fading 
causes initially bright objects to contribute to lower-L 
bins at later ages, thus if the epoch of fading domi- 
nates over the epoch of zero-age luminosity, we would 
expect the observed HII LF to be steeper than the 
unevolved one. We found that for f{t) of the form in 
equation H, and reasonable ranges for Wup, ims('Tiup)j 
and r], that essentially no change in HII LF slope is 
expected. There is an exception for initially shallow 
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Fig. 1. — Evolution of ionizing photon emission rate go (sec ^) from Schaerer & de Koter (1997). 
models are, from bottom, 20, 25, 40, 60, 85, and 120 Mq. 
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slopes of < 77~^ + 1, for which the observed slope in- 
deed steepens to a value a = i]^^ + I. This critical 
value of the slope is probably in the range 1.2 - 1.4, 
which is lower than most measured HII LF slopes, 
hence it is likely that this effect is unimportant. These 
effects are described in further detail in Oey & Clarke 
(1997). Thus in the saturated regime, no change in 
HII LF slope is reasonably expected. 



3. Monte Carlo Models 

The statistical uncertainties that dominate the un- 
saturated clusters render an analytic investigation of 
these objects unfeasible. This is due primarily to 
the fact that many different combinations of stellar 
masses can contribute to any given luminosity, and 
there is no simple way to relate the total nebular lu- 
minosity to the stellar composition, and hence, lumi- 
nosity evolution. We have therefore modeled the H II 
LF with Monte Carlo simulations to investigate its 
behavior in the unsaturated regime. 

Our model draws the integer number of ionizing 
stars from a power-law distribution, so that N{N^,) dN^, 
is the number of clusters with quantities of stars in 
the range TV, to N^ + dA^, : 



N{N^) dN^ = N-^' dN: 



(4) 



We consider "clusters" containing a minimum of one 
star, with the probability of larger clusters decreas- 
ing with A^, according to the exponent /3. Our models 
use the RANI portable random number generator de- 
scribed by Press et al. (1986), to generate the random 
deviate x, having values between and 1. The dis- 
tribution in N^, given by equation |j is then obtained 

by 

1/(1-/3) 



iV, 



<l-/5) 



(5) 



where iV, is assigned to the nearest integer. Since the 
saturated end of the H II LF does not change slope, 
as outlined in the previous section, we adopt /3 = 2 as 
a typical value found for HII LFs in nearby galaxies 
(e.g., KEH). Unlike MW97, we do not consider a dis- 
tribution that is truncated at an upper limit in N^, for 
our default models, although this condition is relaxed 
in § 3.3. 

For each cluster, we then draw iV, stars from a 
stellar mass distribution truncated above a mass limit 



n{m) dm oc 



/ m \ 

\ 771up / 



dm 



(6) 



where n{m) dm, is the number of stars in the range 
m, to m + dm. These are generated with an algo- 
rithm analogous to that for N„ . Equation ^ describes 



a power-law for m <C m^p, and is the same ana- 
lytic form used by MW97 for their N^, distribution 
and IMF; we refer to this form as a truncated power- 
law, following their nomenclature. We use a Salpeter 
(1955) power-law slope of 7 = 2.35, and adopt an 
upper-mass limit iriup = 100 Mq. We include stars 
down to a lower-mass limit mio = 17 Mq (see below). 
The Ha luminosity I associated with individual stars 
is dependent on the mass and main-sequence lifetime 
inis, for which we also use power law parameteriza- 
tions: 

mocl^ , (7) 

normalized so that /up = limup) = 3 x lO'^^ergs"^; 
and 

m oc t-^i , (8) 

normalized so that tms("T-up) = 2.8 Myr (Schaerer 
etal. 1993). For stars in the range 17 < m < 100 Mq, 
we find a value of d = 0.7 from Schaerer et al. (1993). 
We caution the reader that the power-law represen- 
tations of equations M and @ are approximations; the 
relations steepen significantly toward lower masses. 
As discussed below, larger values of S imply a greater 
sensitivity of I to m, thereby increasing the range of L 
in the unsaturated regime, for a given range of m. We 
adopt S ^ 1.5, which yields a reasonable log Zio ~ 37.3 
for mio = 17Mq. For illustrative purposes, it is help- 
ful to show a larger dynamic range of unsaturated 
objects, though we caution that in reality the range 
of these objects is likely to be somewhat smaller for 
our mass limits, since the m — I relation flattens to- 
ward higher masses (e.g., Leitherer 1990). We have 
adopted the value of mio to coincide roughly with the 
value at which the I — m relation turns over more 
steeply, since stars below ~ 17 Mq do not contribute 
strongly to the ionizing fiuxes. There is also a rela- 
tive increase in tmsim) at lower masses, but not as 
significant as the change in l{m). It is useful to bear 
in mind that the representative values of 5 and d will 
increase toward lower masses. 

We consider two different formulations for the stel- 
lar ionizing rate qo with time. The first simply ap- 
proximates that qo remains constant until tms and 
is zero thereafter. The second formulation accounts 
for evolution of qo during the stellar main-sequence 
phase. Figure |l| shows qo for stars of different ini- 
tial masses as a function of time, from the models of 
Schaerer & de Koter (1997). These models suggest 
that I remains fairly constant until a time after which 



it decreases roughly along the relation. 



^ — ^up 



t'msV^upJ 



-c 



(9) 



where C ~ 5, and t is the age from birth. Given the 
uncertainties in massive star evolution and ionizing 
fluxes, we adopt the crude approximation of constant 
/ until the time of intersection with equation 0, af- 
ter which we adopt I from that relation. One conse- 
quence of this formulation is that stars with initially 
different luminosities can have the same I during cer- 
tain of their fading periods. Figure |l| suggests that 
this may not be far from reality. We caution that 
the assumption of coeval star formation within indi- 
vidual nebulae is likely to be an oversimpliflcation, 
and that the stellar age spread is typically similar to 
the main-sequence lifetime of the most massive stars 
(e.g., Massey etal. 1995). Therefore, the early neb- 
ular luminosity evolution could potentially increase 
somewhat before the onset of fading. Speciflc model- 
ing of star formation scenarios would be necessary to 
investigate this effect. 

Our adopted lower-mass cutoff introduces artifi- 
cial fluctuations that can be seen in all our models 
at low L. The quantized nature of the stellar ioniz- 
ing sources, combined with the truncation of stellar 
masses below ttiio, causes artificial kinks in the dis- 
tribution of nebular luminosities at low L. For ex- 
ample, the most prominent of these features appears 
at logL = 37.6 in our models, or 0.3 dex above the 
lower-L cutoff of log/io = 37.3. This dip results from 
the fact that all nebulae at L ;^ 37.6 contain single 
stars, selected according to the IMF, whereas two-star 
nebulae begin to contribute above this luminosity. To 
more clearly demonstrate the features in the H II LF, 
we therefore include 10'* clusters in our simulations, 
although this number is higher than is seen in most 
nearby galaxies. Finally, we truncate our model HII 
LFs below the minimum luminosity l\o, which is that 
due to a single star of mass m\o. Although objects 
originating at higher L evolve to luminosities below 
ho for cases incorporating main-sequence evolution, 
the omission of objects originating below ^lo renders 
the model populations incomplete below that value. 

3.1. Single Burst 

Figure shows two models for the evolution of an 
H II region population created in a single burst, with 
6 = 1.5, d = 0.7, and 17 < m < 100 M©. The left 



column of Figure || (panels a — c) shows model Al, 
with constant I during the main-sequence; while the 
right column (panels d — f) shows model A2, with l{t) 
as described above, to include main-sequence fading. 
The zero-age H II LFs are given in the top panels, with 
the middle and bottom panels showing ages of 4 and 7 
Myr, respectively. Note that the two zero-age models 
(top panels) are statistically identical. To aid com- 
parison of the models, the dotted line corresponding 
to Lsat = 3 X 10"^^ ergs^^ is shown for reference. 

As found in the simulations of MW97, the zero- 
age models in Figure (top panels) show a flattening 
of the HII LF slope below a saturation luminosity. 
For the adopted stellar parameters, this transition 
from saturated to unsaturated populations occurs at 
logLsat = logZup = 38.5. Equation shows that the 
extent of this unsaturated regime is determined by 
S: a larger S implies a greater sensitivity of I to m, 
thereby increasing the range in L of this population, 
for a given range of m. This effect also causes the 
slope in the unsaturated regime (L < Lgat) to be de- 
pendent on 6 because the scatter in L contributed by 
individual stars is larger for large S, and this scat- 
ter determines the flattening in slope. Hence, larger 
values of S will induce a flatter slope. As mentioned 
above, our adopted value of S = 1.5 is an overestimate 
for the more massive stars in the IMF, hence the un- 
saturated slope should be somewhat steeper near Lgat 
than is seen in our models. 

H II regions with L > Lsat are not sensitive to the 
m — I relation since statistical effects are unimportant 
in this regime. The slope of the HII LF in the sat- 
urated regime is (3, the same as that of the parent 
distribution in N^ (equation ^, as described in §2. 

Figures |p and c show Model Al at ages of 4 and 7 
Myr respectively. These ages should not be regarded 
too literally, since the evolutionary state at a given 
age is dependent on equation ||, which is an approxi- 
mate parameterization, as well as equations and R. 
For reference, tms('7iio) = 9.7 Myr for our adopted 
d — 0.7. It is immediately apparent that the inflec- 
tion in slope moves toward progressively lower lumi- 
nosities with age. This is primarily a reflection of the 
group evolution of the saturated population to lower 
luminosities; the entire distribution at L > Lsat in 
Figure ^a is moving toward lower luminosities in Fig- 
ures |p and c. For cases such as this, in which the 
zero-age Lgat = ^up, the slope inflection at later ages 
now corresponds to I of the most massive star still 
present in these objects. More generally, the rate at 



which the saturated objects fade in L has a depen- 
dence of t^^'^, as determined by equations and ||. 

Model A2 includes main-sequence evolution of the 
ionizing fluxes. We caution that the models includ- 
ing the main-sequence evolution should be interpreted 
qualitatively, owing to the crude parameterization of 
the stellar evolution. The principal feature is that 
Model A2 clearly evolves much faster than Model Al. 
The other important effect is that noticeable, sharp 
structure in the HII LF becomes apparent near the 
transition from saturated to unsaturated objects, as 
seen in Figure He as contrasted to Figure pp. This re- 
sults from the fact that clusters dominated by a few, 
high-mass stars are now capable of evolving to lower 
L than clusters dominated by many, unevolved, low- 
mass stars. In other words, stellar evolution allows 
the low-L end of the saturated population to over- 
take, toward lower luminosities, the high-L end of the 
unsaturated population. The details of this structure 
in the HII LF are dependent on the details of the 
main-sequence luminosity evolution. 

3.2. Continuous Creation 

Figure shows simulations for a population where 
the H II regions are being continuously created (Model 
Bl). This model is the same as that in the sin- 
gle burst case with no main-sequence luminosity evo- 
lution (Model Al), but now with each new cluster 
formed at a uniform incremental time after its prede- 
cessor. The total distribution of initial luminosities, 
logLoj of all the objects, which represents the parent, 
unevolved HII LF, is the same as the zero-age dis- 
tributions in Models Al and A2 (Figures |p and d). 
Figure ||a shows the distribution at an age of 6 Myr 
after the onset of continuous creation, but before the 
population has achieved its steady-state form. The 
steady-state HII LF, which occurs at i > ims("T'io)) is 
shown in the Figure |3p. 

We may consider the HII LF for continous cre- 
ation to be a sum of individual bursts with a contin- 
uous age distribution. Since the single burst HII LF 
evolves predominantly as a progression of the satu- 
rated objects, with slope /?, to the lower luminosities 
of the unsaturated regime, we therefore expect that 
the slope in the region L < Lgat will be intermediate 
between its flatter value in the unevolved case and 
/?. Figure |^6 shows that this is indeed the case. We 
can also see in Figure 0a, that the intermediate evolu- 
tionary stage shows the region of steepening moving 
in from higher luminosities, since the low-luminosity 
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Fig. 2. — Model HII LFs for a single burst. Model Al (panels a — c) considers constant / over main sequence 
lifetimes; Model A2 (panels d — f) considers main sequence fading of I as described in §3. The dotted line marks 
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bins are not yet affected by saturated, evolved ob- 
jects. We also generated a Model B2, incorporat- 
ing main sequence luminosity evolution, analogous to 
Model A2. The steady-state form of that model also 
shows an intermediate- value slope at L < Lgat, but 
slightly flatter than its analog in Model Bl because 
of the structure present at i < Lgat in the burst case 
(Model A2). 

3.3. Unsaturated Population 

In Figure 0, we show models for a population of 
only unsaturated clusters. These are drawn from the 
same power law distribution in N^ (equation |j) as the 
previous models, but only for N^ < 10. For our de- 
fault stellar parameters with 17 < m < 100 M©, the 
saturation transition in A^, is roughly between 10 and 
20 stars. We use the same stellar parameters as in 
the previous models, now with a total of 3000 clus- 
ters. The left side of Figure ^ shows the simulation 
for the unevolved population (a) , and continuous cre- 
ation (6); while the right side shows the single burst 
case at ages 4 Myr (c) and 7 Myr (d). These models 
do not include main sequence luminosity evolution, 
whose consequences are analogous to those discussed 
for Models A2 and B2 above. 

It is apparent in Figure H that the primary dif- 
ference in the HII LF between this population and 
the complete power- law sampling in A^,, is that the 
high-L tail of the unsaturated objects falls off very 



steeply. The highest-luminosity, unsaturated clusters 
result from a chance assemblage containing a dispro- 
portionate number of high-mass stars. Therefore, the 
distribution of these objects must fall more steeply 
than that of a saturated population, owing largely to 
the scarcity of the high-mass stars. The actual slope 
in this regime is once again determined by the detailed 
statistics of the cluster compositions, and equations^ 
and 0. 

At L < Lsat, the behavior of these unsaturated 
populations is identical to that seen in the models 
with a full distribution in A^*, since, as discussed 
above, this range in L is dominated by the unsatu- 
rated objects. It is only at the high-i limit itself, in 
the evolved case of Figure ^, that the sudden drop 
due to sheer lack of objects, effectively generates a 
steep slope. Note that at luminosities below this high- 
L drop-off, the H II LF in this model recovers a slope 
of /?, analogous to Figure gc. 

4. Observed HII LFs 

The features and behavior found in the simulated 
HII LFs can provide a useful interpretation of the 
patterns found in observed HII LFs. 

4.1. Arm and Interarm HII LFs 

In Figure o, we reproduce the published H II LFs 
for arm and interarm populations in six grand de- 
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Fig. 4. — Models with the same parameters as in Figures and 0, but for populations with N^ < 10. (a) Unevolved 
population; (6) steady-state, continuous creation; (c) single burst at age 4 Myr; (d) single burst at age 7 Myr. To 
aid comparison, these models are plotted on the same scale as Figures |^ and |[ 



sign spirals. The solid histograms show the arm HII 
LFs, and the dot-dashed lines show the interarm HII 
LFs. In each case, the peak of the interarm HII LF 
occurs at a fainter luminosity than that for the arm 
regions, suggesting that the interarm H II LF mirrors 
that of the arms at lower L. This shift in peak oc- 
curs despite the smaller total numbers of H II regions 
in the interarm areas. Owing to incompleteness ef- 
fects, we might otherwise expect the smaller numbers 
to induce a turnover at, if anything, higher L. We 
do caution that blending of nebulae (Deharveng et 
al. 1988; KEH) in the more crowded arm populations 
could also introduce a flattening in the observed H II 
LF at lower L. However, as pointed out previously 
(Rand 1992; Knapen etal. 1993; KEH), the highest- 
L objects in the interarm regions are apparently al- 
ways fainter than those in the arms. This is again 
consistent with a group shift in luminosity. 

These features exhibit the behavior of our model 
HII LFs for the burst scenario. This is consistent 
with the expectation that the arm nebulae are largely 
coeval, and are undergoing a current burst of cluster 
formation; whereas the interarm regions represent, on 
average, an older population remaining from previous 
star formation. As seen in Figure 0, the inflection in 
the H II LF occurs at progressively lower values with 
age. Note that even if incompleteness affects the lo- 
cation of the peak for the arm objects, which we take 
to reflect the unevolved population, the relative shift 
in the peak for the interarm objects is still consistent 
with evolution, owing to the steeper slope at lower L. 
This is because the low-L bins are more robust to in- 
completeness for a steep slope than for a shallow one, 
since the proportion of missing objects is the same in 
each case. Therefore, while the absolute location, or 
existence, of the H II LF peak may be suspect due to 
observational limitations, the observed shift in peak 
for the interarm regions is fully consistent with our 
model because a steeper slope is expected at low L 
(cf. Figure gd and /). Indeed, given that the HII 
LF peak is often substantially lower in log A^ (log L) 
for the interarm regions, this suggests that the actual 
peak of the distribution has moved to luminosities 
much lower than the empirical completeness limit, as 
is illustrated, for example, in Figure g/. 

The aging effect can also explain the steeper HII 
LF slopes that are sometimes reported for the inter- 
arm regions as compared to the arm regions of spiral 
galaxies (Point 2 in §1). If the actual nebular lumi- 
nosity function conforms to the two-slope shape seen 



in the models, and if the two slopes are not fitted 
separately, then the net effect is that the measured 
slope will be intermediate between the two. In our 
single burst model, the location of the break moves 
to lower L. Therefore, over a fixed range of L, the 
contribution of the shallower slope to the composite, 
measured HII LF slope will be less for an evolved 
burst population, than for an unevolved one. As is 
apparent in Figure 0, it is difficult to pinpoint the lo- 
cation of the slope change. We mention once again 
that our simulations show a flatter slope in the un- 
saturated regime than would be expected, owing to 
the parameterization of equation 0. Furthermore, the 
observed location of Lsat will also vary with galactic 
internal extinction, which most HII LFs, including 
the ones in Figure || do not take into account. With 
these difficulties, we suggest that the slopes of the in- 
terarm regions do appear to be the same as those in 
the saturated regime of the arm populations. 

Thus the differences observed in the HII LFs for 
arm and interarm regions are not necessarily due to 
changes in star formation environment, but are likely 
to be affected by evolution. The time scale between 
the passage of spiral density waves is of order ~ 40 
Myr (e.g.. Rand 1993) at any given location in the 
disk, which is much longer than the observable life 
expectancy of the H II regions. Thus, the spatial and 
luminosity distribution of nebulae between the arms 
does imply that there is at least a certain amount of 
bona fide interarm star formation taking place. How- 
ever, based on the arguments presented in §1, we feel 
that an aging population is likely to be an important 
component in the interarm regions, and that the be- 
havior of the H II LF is consistent with this interpre- 
tation. With the present data and the observational 
uncertainties mentioned above, combined with uncer- 
tainties in predicted go a-nd its evolution, we have not 
attempted to place meaningful age constraints on the 
evolved population. However, it will be interesting to 
re-examine the data more closely, and search for an 
evolutionary progression in the H II LF with distance 
behind the leading edge of the spiral arms. We hope 
to carry out such an investigation in the near future. 

4.2. The Hubble Sequence 

The correlation of earlier type galaxies with steeper 
slopes (Point 3 in §1) can be explained by systemati- 
cally lower cutoffs in A^* in those clusters. KEH have 
shown compelling evidence that these galaxies have 
fewer stars per cluster, since the mean nebular lu- 
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minosities in the early type spirals are substantially 
lower than those in late types, despite the larger total 
numbers of H II regions in the earlier types. There- 
fore, unless the stellar ionizing population varies dras- 
tically, it is probable that the size of the parent clus- 
ters must themselves be smaller on average in the 
earlier type galaxies. The same conclusion has been 
found for Sa galaxies (Caldwell etal. 1991), although 
the total numbers of H II regions are again smaller in 
those galaxies. 

The decrease in typical N,, per cluster with Hubble 
type (KEH) could easily result from a decrease in the 
upper cutoff of the N^ power-law. The earliest type 
galaxies would have populations composed only of un- 
saturated objects, as would correspond to Sa galax- 
ies, where the maximum log L is typically in the range 
38.0 - 38.5 (Caldweh etal. 1991). This is the same 
range in values as the estimates for ^up, as mentioned 
earlier, hence it is indeed the regime of unsaturated 
objects. As seen in §3.3, the models that describe this 
population show a steeper high-luminosity tail above 
isat, than do the populations with saturated objects. 
Indeed, the slopes in Sa galaxies have been found to 
be even steeper than those for Sb - Sc types: the 
typical HII LF slope a ^ 2.6 for Sa galaxies (Cald- 
well etal. 1991), as compared to a ~ 2.0 for Sb - Sc 
galaxies (KEH; Banfi etal. 1993). The model in Fig- 
ure^, for continuous creation of unsaturated objects, 
would therefore be the most appropriate of our sim- 
ulations to describe the Sa HII LFs. The form of the 
Sa H II LFs shown in Figure 5 of Caldwell et al. (1991) 
show a consistent resemblance with our models in Fig- 
ure 0. Thus, ironically, although these nebulae are 
dominated by unsaturated objects whose characteris- 
tic slope is flatter than the saturated case, the high-L 
drop-off will induce a steeper measured slope in these 
populations. 

The Sb - Sc galaxies would then show a slope inter- 
mediate between the Sa and Im types (cf. Figures 46 
and 36) , as the upper cutoff in N^ approaches the ex- 
clusively unsaturated regime. Indeed, KEH suggest 
that the HII LF in Sb - Sc galaxies is truncated at 
an upper-i limit. This is quantitatively supported by 
MW97, who also find evidence for a truncated H II LF 
in the Milky Way and Sb-Sc galaxies, but not in ir- 
regular galaxies. Since this limit is still at L above the 
regime of unsaturated objects, these findings are con- 
sistent with our interpretation of a progressively lower 
upper- iV» cutoff with earlier Hubble type. There may 
be exceptions caused by temporarily enhanced star 



formation, in which large numbers of saturated clus- 
ters are formed. This scenario can explain the larger 
range in HII LF slope seen in earlier type galaxies 
(KEH). Indeed, a slope steepening at high L reported 
in some earlier type galaxies could well correspond to 
a cutoff; indeed, KEH note that it is often impossible 
to distinguish such a turnover from a cutoff. Thus the 
steepening in slope for disk galaxies across the Hubble 
sequence can be explained by a simple trend in upper 
cutoff of the iV* distribution, and need not invoke a 
change in the parent slope /3 of A^, . 

4.3. The Slope Transition 

In §3, we denoted the luminosity of the slope break 
in the initial HII LF by Lgat and noted that, as in the 
models of MW97, this often corresponds to the lu- 
minosity contributed by most luminous star in the 
IMF, /up. We however caution that this correspon- 
dence is dependent on the stellar parameters. As 
an example. Figure g exhibits the HII LF resulting 
from two models with identical IMF slope and m — I 
relation, but with different stellar mass limits. Fig- 
ure |a, with 3 < m < IOOM0, exhibits an HII LF 
that more closely resembles in form the simulation 
by MW97. There are two principal slope breaks ap- 
parent: one near l^p] and one near I, below which 
the slope appears quite flat. The existence of a min- 
imum cluster membership number, N^ — 1, causes 
this flattening, or in some cases even positive slope, 
for luminosities below the expectation value of L con- 
taining single stars, i.e., the average stellar luminosity 
L At I < L < l^p, there is an intermediate slope in 
the HII LF that joins the flat regime to the satu- 
rated regime of slope /3. A small transition region of 
<^ 0.3 dex is noticeable in this intermediate-slope re- 
gion near both / and /up, which results from the form 
of the IMF adopted in equation ^. With instead a 
sharp IMF truncation at niup, the slope breaks would 
occur cleanly at I and l^p- 

One can now see that certain stellar mass ranges 
and values of S can change the relative location of 
the two breaks such that one or the other could dis- 
appear. The model in Figure 06, with 17 < to < 
300 M0 , shows an example where the l^p slope break 
is virtually undetectable at logZup = 39.2. The 
smaller difference between I and ^up has compressed 
the intermediate-slope range in L-space, thereby steep- 
ening its slope to a value closer to (3, which represents 
the maximum slope possible in this regime. An ob- 
served H II LF similar to Figure pp would most likely 
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be measured to have only one slope in the regime 
L > T. Note that the models in Figure ^ are prob- 
ably unrealistic, since 5 = 1.5 is not representative 
over such large mass ranges, and we only show these 
for illustrative purposes. However, these examples 
demonstrate that breaks in observed H II LFs should 
be interpreted with caution, and may only provide a 
lower limit to the most luminous star in the IMF. Fur- 
thermore, the tendency of the highest-mass stars to 
approach similar ionizing emission rates go (e.g., Lei- 
therer 1990), compounds the insensitivity of ^up to 
'Tiup- We therefore conclude that, unfortunately, the 
HII LF is not a useful probe of TOup. 

However, the insensitivity of Zup to TO,up suggests 
that the high-L end of the H II LF should be a reliable 
indicator of /3 for saturated nebular populations, since 
it diminishes the likelihood of a varying TOup as an 
explanation for the observed variations in the H II LF. 
The principal caveat is that a truncation at maximum 
iV, not be confused with a steep slope of /3, as may 
have been done for Sa galaxies in the past (see above) . 
It is also apparent that, since the form of the H II LF is 
sensitive to the stellar parameters, that demonstrable 
variations in the unsaturated regime of the zero-age 
or continuous creation H II LF could possibly indicate 
substantial variations in the stellar parameters of the 
hottest stars. 

We note that stochastic variations in TTiup within a 
galaxy will not change the location of the break if the 
distribution of TO,up is gaussian. We confirmed this 



by drawing rriup from a parent normal distribution of 
cr = 15 Mq and the same mean value of 100 Mq. The 
resulting H II LF is indistinguishable from that with 
a constant TTiup- 

Empirically, a flattening at lower L is often seen 
in many galaxies (Point 1 in §1), but at varying crit- 
ical luminosities. While the models for a single burst 
show that the transition point is not necessarily fixed, 
Lsat cannot be greater than Zup- At present, the value 
of ^up, the Ha luminosity associated with TTiup, is not 
well-determined. Our models have adopted a value 
of Zup = 3 X 10"^^ ergs^^, which is uncertain within 
a factor of a few, and is a high estimate. Most pub- 
lished HII LFs are not corrected for internal extinc- 
tion within the galaxies, which will affect the location 
of the break. As can be seen in Figure 86 of Walterbos 
& Braun (1992), extinction correction for the highly- 
inclined case of M31 can shift the location of the ap- 
parent turnover by an order of magnitude. It is also 
important to distinguish between a real slope change 
and artificial slope flattening induced by limitations 
such as incompleteness, and the more insidious prob- 
lem of object blending, mentioned above, for more 
distant galaxies. In view of the difficulty in pinpoint- 
ing the location of a slope change, and the variations 
induced by reddening and possible aging effects, we 
suggest that the published H II LF slopes are thus far 
consistent with our models. 

We note that KEH modeled the lower-L end of 
the H II LF by combining the m — I relation with the 
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IMF for single-star H II regions. This representa- 
tion is somewhat different from ours, since we include 
multiple-star, low-L objects that will contribute to 
the same luminosity bins occupied by bright, single- 
star nebulae. Our slope in this regime will there- 
fore tend to be slightly flatter (more positive) than 
theirs. This will affect estimates for the fraction of 
total galactic ionizing emission due to faint H II re- 
gions below empirical completeness limits, such that 
those using the KEH model are an upper limit for 
similar stellar parameters. 

Additional models of the HII LF by Feinstein 
(1997) came to our attention during the refereeing 
process. These models consider only nebular popula- 
tions with luminosities scaled from objects with a sat- 
urated IMF, but with mio extending down to SM©. 
We emphasize that, since unsaturated clusters are not 
considered, the slope turnover seen in those models is 
therefore not caused by the same mechanism seen in 
our work, or that of KEH and MW97. Rather, the 
slope break seen in Feinstein's models result from the 
adopted lower cutoff in the L distribution, combined 
with an extended creation period. The zero-age, sat- 
urated H II LF is described by the power- law /3, which 
simply shifts to lower luminosity with advancing age, 
as exhibited by the saturated tail of our models in 
Figure g. Feinstein's continuous creation models are 
therefore a superposition of many power-laws with a 
continuous range of ages. The lower-L cutoff seen in 
those models therefore induces a flattening in the H II 
LF below that luminosity. This is best illustrated in 
the right panel of Feinstein's Figure 3, which varies 
the lower-L cutoff according to assumed total stellar 
mass, for continuous creation. The slope break oc- 
curs at the lower-i cutoff, hence, as pointed out by 
the author, this parameter determines the location 
of the break. We caution, however, that the lower-L 
cutoffs used in Feinstein's models are necessarily ar- 
tificial, owing to the consideration of only saturated 
H II regions. Therefore, the slope breaks seen in these 
models must be regarded as artificial to the same ex- 
tent. If, however, nebular creation scenarios exist 
with such a strong lower-L cutoff, Feinstein's mod- 
els demonstrate another mechanism which can also 
induce a slope break. 

Previous work has sometimes suggested an empiri- 
cal slope change at high luminosities of log L ^ 38.6 or 
even higher (Rozas etal. 1996a; Rand 1992; KEH). If 
a break at these high luminosities is real, it is unlikely 
to result from the stochastic effect below Lgat- How- 



ever, if the high-L slope is steep enough, it could well 
be a manifestation of the maximum TV, cutoff in those 
galaxies, which is expected to vary in the scenario we 
presented in §4.2. It could also result from a change in 
Lsat, although this appears unlikely, as argued above. 
But factors other than those that we have considered 
can affect the form of the H II LF as well. For exam- 
ple, Rozas etal. (1996a, b) suggest that such a break 
at high L results from a transition between radiation- 
bounded nebulae at fainter luminosities, to density- 
bounded objects at brighter luminosities. Another 
potentially important effect is a differential extinction 
pattern among the nebular population, and of course, 
the usual differences in IMF slope and star forming 
conditions. Given the currently available data, the 
systematic presence of an additional, high-luminosity 
slope change at a fixed L remains controversial. 

5. Conclusion 

We have modeled the HII LF with Monte Carlo 
simulations that draw A^, for each cluster from a 
power-law distribution of constant slope, and for each 
cluster, m from a Salpeter IMF truncated at TTiup = 
100 Mq. As found by MW97, this produces shallower 
slopes corresponding to the population dominated by 
unsaturated objects, and a steeper, high-L side dom- 
inated by saturated ones. The saturated slope is /3, 
the slope of the A^, distribution (equation Q). The 
transition luminosity Lgat in the unevolved HII LF 
may correspond to /up, the Ha luminosity due to a 
single star of mass rriup, but this depends on the as- 
sumed stellar parameters. The observed HII LFs do 
often show flatter slopes at lower luminosities that are 
consistent with the transition at Lgat (e.g., Walterbos 
& Braun 1992; Scowen, Dufour, & Hester 1992). 

We have hitherto not discussed the form of the H II 
LF below l\o, which is unexplored in our models. Since 
the slope S of the m — I relation (equation |^) steep- 
ens, the HII LF will continue to grow flatter for the 
same A^, power law, extending to lower luminosities. 
At logL <^ 36, contamination from planetary nebulae 
may become important (Walterbos & Braun 1992), 
while confusion and background diffuse Ha emission 
will undoubtedly limit the HII LF in this regime. 

Our simulations show that for a single-burst case, 
the slope transition moves to lower L as the popula- 
tion ages and the saturated objects grow fainter (Fig- 
ure ||) . Hence, the steady-state model for continuous 
creation shows a slope in the region L < Lgat that is 
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intermediate between its value in the unevolved HII 
LF and the slope (3 of the saturated objects (Figure ||). 
Models that include main-sequence evolution for stel- 
lar ionizing luminosities show a faster evolution than 
those that do not, with a noticeable peak around the 
transition between saturated and unsaturated popu- 
lations (Figure He). The steady-state model for con- 
tinuous creation incorporating stellar evolution yields 
a slightly shallower slope than that for the case with- 
out stellar evolution, but still steeper than the un- 
evolved case. Models for a population of only unsat- 
urated clusters with N^, < 10 show a much steeper 
drop-off at L > Lsat (Figure H), where the HII LF 
is normally dominated by saturated objects. The be- 
havior at L < Lsat is generally the same as that for 
the models with fully populated N^,, with the excep- 
tion of a high-L drop-off in an evolving burst. 

The patterns seen in observed HII LFs in nearby 
galaxies can be explained by the behavior of these 
models. If H II regions in spiral arms of galaxies repre- 
sent a current burst of coeval star formation, then the 
steeper slopes seen in interarm vs. arm regions can 
be explained as an aging burst compared to a current 
burst population. Since single slope measurements 
will yield a composite of the saturated and unsatu- 
rated slopes, the zero-age H II LF will show a flatter 
slope owing to the larger range of unsaturated ob- 
jects. The lower peak L and maximum L seen in the 
interarm regions are also consistent with a dominant 
contribution from an aging population. 

The steeper slopes seen in earlier galactic Hubble 
types can be explained by the lower cutoff in max- 
imum iV,. As seen in the models for populations 
of only unsaturated clusters (Figure 0), the high-L 
tail is much steeper than in populations with higher 
A^* . A truncation at the high end of the H II LF has 
been suggested by KEH and MW97. Since the mean 
A^* does apparently decrease in earlier type galaxies 
to the values in unsaturated clusters (KEH; Caldwell 
etal. 1991), the steeper HII LF slopes are likely to 
simply result from a lower maximum A^*. A slope 
dropoff therefore may also be seen at L correspond- 
ing to the maximum A^*. 

Thus, we suggest that the current data are consis- 
tent with evolutionary effects and the upper- A^* cutoff 
being the primary determinants of the shape of the 
HII LF. Indeed, it may be unnecessary to invoke any 
other parameters for the nebular and cluster prop- 
erties to explain the behavior of the HII LF as de- 
scribed in §1. If this is true, it is rather remarkable, 



for example, that grand design spirals with strong 
star formation induced by tidal interactions, would 
show the same power-law distribution in N^, as those 
found in normal disk and irregular galaxies. However, 
we again acknowledge that other, unexplored effects, 
such as gas density bounding, differential reddening, 
and variations in IMF and other stellar parameters, 
can influence the form of the H II LF. Nor can we rule 
out suggested environmental influence like variations 
in slope of the molecular cloud mass spectrum. But 
thus far, it appears that the power-law slope f3 of A^* 
might vary even less than was suggested by the ob- 
served slopes of the H II LF, even though the rate and 
total star formation can vary by one to two orders of 
magnitude. Further detailed observations of the HII 
LF, and re-examination of those in the literature, are 
of great interest in evaluating the possibility of a uni- 
versal /3 ~ 2. Suggestively, Caldwell etal. (1991) find 
that the luminosity function of knots seen in the B 
continuum of four of their Sa galaxies do show such a 
power-law exponent ^ 2. 

As is done with the luminosity functions of distant 
galaxies, it may thus be possible to use the form of the 
H II LF to understand features in the most recent star 
formation history (< 10 Myr) of their host galaxies, 
since individual bursts can leave their imprint on the 
HII LF. Our models, for example, show that an aged 
burst population is distinguishable from continuous 
creation of low-A^, clusters (cf. Figures Ec,/ and Ep). 
Feinstein (1997) also demonstrates effects of different 
creation histories. At present, it is difhcult to make 
useful quantitative constraints, owing to present con- 
straints on stellar and cluster parameters, but it is 
at least possible gain qualitative insight on the recent 
star formation history. 

We thank Rob Kennicutt, Johan Knapen, and 
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manuscript. We also thank Johan Knapen for access 
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